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Abstract 

We present detailed calculations of the quasinormal modes of Reissner- 
Nordstrom black holes. While the first few, slowly damped, modes depend 
on the charge of the black hole in a relatively simple way, we find that the 
rapidly damped modes show several peculiar features. The higher modes gen- 
erally spiral into the value for the extreme black hole as the charge increases. 
We also discuss the possible existence of a purely imaginary mode for the 
Schwarzschild black hole: Our data suggest that there is a quasinormal mode 
that limits to tuM = — 2i as Q —* 0. 
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I. PERTURBING THE REISSNER-NORDSTROM BLACK HOLE 



Quasinormal modes of black holes have been studied ever since the seminal work of 
Vishveshwara [[J and Press 0] in the early 1970s. It soon became evident that exponentially 
damped mode-oscillations will dominate most processes involving perturbed black holes (see 
H for references). This means that the quasinormal modes provide a unique opportunity to 
identify a black hole, a possibility that hopefully will become reality when large-scale laser- 
interferometric detectors for gravitational waves come into operation in the near future 
In order to extract as much information as possible from a gravitational-wave signal it is 
important that we understand exactly how the quasinormal modes depend on the parameters 
of the black hole. 

The parameters of main astrophysical importance are the black holes mass and angular 
momentum. That is, the Kerr solution is the most relevant one from an astrophysical 
point of view. The solution that describes an electrically charged, nonrotating black hole 
- due to Reissner and Nordstrom — is of less direct importance because it seems unlikely 
that black holes with a considerable charge will exist in the Universe. Nevertheless, the 
Reissner-Nordstrom solution has several interesting features that warrant a closer inspection. 
The most intriguing one concerns the possible conversion of electromagnetic energy into 
gravitational energy and vice versa: In a charged environment an electromagnetic wave will 
inevitably give rise to gravitational waves. The Reissner-Nordstrom metric provides the 
simplest framework for studies of this effect. 

For this and several other reasons there has been a number of studies of perturbed 
Reissner-Nordstrom black holes. The equations governing a weak (massless) field in the 
geometry of an electrically charged black hole were first derived by Zerilli || and Moncrief 
||. Basically, these equations can be i) split into axial and polar perturbations (also known 
as odd and even parity, respectively) and then ii) reduced to two decoupled wave equations 
in each case. The final wave equations describe two variables \&i and ^2, from which all 
components of the electromagnetic field and the perturbed metric can be reconstructed. 
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In general, both and \&2 correspond to a combination of electromagnetic waves and 
gravitational ones, but in the limiting case of an uncharged black hole the two functions 
reduce to pure electromagnetic and gravitational waves, respectively. 

Originally, the perturbation equations were used to study the stability of the metric . 
Later the equations were used to investigate the already mentioned conversion of electro- 
magnetic energy into gravitational [0-f| . There has also been several studies of quasinormal 
modes for Reissner-Nordstrom black holes ||lO|-[r^| . But although the methods used in those 



studies provide accurate numerical results there is still need for more information. Specif- 
ically, one would like to know what happens to the quasinormal modes as the black hole 
becomes extremely charged. There are two parts to this problem, that seemingly must be 



studied separately. In the typical case each of the two horizons, at r± = M ± \JM 2 — Q 2 
(where M is the mass and Q < M the charge of the hole), correspond to a second order 
singularity in the linearized equations. But in the extreme case these two singularities merge 
(at r = M) into a single one of fourth order. This means that the methods that have been 
devised to determine the entire spectrum of quasinormal modes break down in the extreme 
limit 



Even though the extremely charged case is only of theoretical interest it is worth con- 
sidering. One reason (apart from natural curiosity) is the fact that the existence of two 
coalescing horizons is a prominent feature also of rapidly rotating black holes. Thus one 
may hope that a study of the mode-behaviour as Q — > M can lead to some insights also 
for the rotating case. Or at least that methods that prove reliable for high charge can be 
adapted to the Kerr case. Recently Leaver's continued fraction method was amended in 



such a way that it could be applied to the case of extreme charge fll3f . Although it seems 



plausible that other methods, such as the numerical integration of Andersson |12|], can be 
extended in a similar way there has been no attempts to do this. 

Anyway, at present the main question concerns the reliability of the various approaches 
for nearly extreme black holes. It is clear that the methods in |TT| , [r2| ] will break down for 



Q = M, and also that the method of Onozawa et al. |T3 can be used only for the extreme 



case. Thus, we do not yet know to what extent the available results for high charge (and 
perhaps also rapid rotation in the case of Kerr [[OJ) can be trusted. The work presented in 
this paper was motivated by a desire to obtain a better understanding of this issue. We also 
wanted to unveil the detailed behaviour of the quasinormal modes as the charge of the black 
hole was increased. Specifically, previous evidence (see Figure 1 in [IT] and Figures 2 and 



3 in |12|]) indicate that this behaviour is somewhat peculiar for the higher overtones of the 
black hole. As we will show in the following section, the behaviour of the highly damped 
Reissner-Nordstrom quasinormal modes is, indeed, very strange. 

II. PRESENTING THE RESULTS 

A. The numerical work 

Since the equations that describe a perturbed Reissner-Nordstrom black hole are available 
in the literature (most notably in Chandrasekhar's exhaustive book [IB]) we will not list them 



here. For the present discussion it is sufficient to know that the equations take the general 
form 

+ [u; 2 -V(r)]* = O, (1) 

where we have assumed that the time-dependence of the perturbation is e~ iuJt . The tortoise 
coordinate, that is defined by 

dr* ( r r 2 ) dr ' ^ 

has the effect that the event horizon of the black hole is "pushed away" to r* = — oo. 

The function *f> can be either an axial (odd parity) or polar (even parity) perturbation 
of the black hole. Also, for each case there are two different functions, \&i and that 
correspond to pure electromagnetic and gravitational waves in the Schwarzschild limit. In 
total there are thus four equations, all with slightly different effective potentials. For obvious 
reasons this difference in the potentials leads to somewhat different quasinormal modes for 



and ^2- Analogously, one might expect the quasinormal modes for axial and polar 



perturbations to be different. But, as pointed out by Chandrasekhar |16[, the mathematical 
theory of black holes is an intricately entangled web where many quantities are related 
(often in a surprising way). Thus it turns out that the quasinormal modes for axial and 
polar perturbations are identical. The underlying reason is that two effective potentials 
can, although different, contain much the same physical information. This means that it 
is sufficient to restrict a study of the Reissner-Nordstrom problem to either the axial or 
the polar case. Once and \&2 are found for one kind of perturbation the corresponding 
solutions for the other case are easily generated . 

The quasinormal modes of a black hole correspond to solutions to Eq. (P that satisfy 
the causal condition that no information should leak out through the event horizon of the 
black hole, and at the same time correspond to purely outgoing waves at spatial infinity. 
This means that a typical mode-solution will behave like 

!e lujr * , as +oo , 

(3) 
e~ tujr * , as — > — oo . 

Since one would expect the black hole to be stable against a small perturbation the mode- 
frequencies should be complex. This implies that an identification of quasinormal modes 
is non-trivial. In order for a mode-solution to be damped with time at a fixed r* the 
frequency u must have a negative imaginary part. But then the corresponding solution 
will diverge both at the event horizon and spatial infinity [cf. Eq. (0)] at a fixed time. 
Several methods have been devised to deal with this difficulty. Most notable are Leaver's 
continued fraction approach |L5| and Andersson's complex-coordinate integration method 



| 17| . These two methods provide highly accurate numerical results also for rapidly damped 
modes. Furthermore, both methods have been used to study the first few quasinormal modes 
for Reissner-Nordstrom black holes. 

As already mentioned in the previous section, it is clear that all existent methods will 
fail for a nearly extreme black hole. The reason is that the two horizons of the black hole 
merge as Q — > M, and this changes the singularity structure of the problem (this is easy to 



see if Eq. ([!]) is expanded in such a way that all derivatives are taken with respect to r rather 
than r*). Thus, the extreme case must be considered separately. This was recently done by 



Onozawa et al. As a test of the reliability of that study, and also to provide a better 
understanding of the behaviour of the quasinormal modes for highly charged black holes we 
decided to make an exhaustive study of the problem. We made detailed calculations [using 
both the continued fraction method fll] and numerical integration [12] to make sure that 
the results were reliable] for the first nine dipole (£ = 1) modes of and the first nine 
quadrupole (£ = 2) modes of ^2- These are the lowest radiating multipoles for each case. 
In general one finds, just as for Schwarzschild black holes ||, that the mode-behaviour is 
similar for higher values of £. The results of this investigation are shown in Figs. [I] and |2|. 
We now proceed to discuss them in more detail. 



B. Slowly damped modes 



For the first few modes the behaviour of the mode-frequencies is readily described. The 
damping rate typically reaches a maximum for Q/M fa 0.7—0.8, and the oscillation frequency 
generally increases with Q. This is clear from the results for n = — 3 in Fig. [I]. Moreover, 
this behaviour agrees with the understanding from the WKB approximation ||10|| . It is 
relatively straightforward to verify that, when the lowest order of approximation is used, 
the WKB formulae suggest the following approximate behaviour (for the slowest damped 
mode) 
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in the limit £ 3> 1. Here we have defined r as the position where the black-hole potential 
attains its maximum value. This means that 2ro ~ 3M + a/9M 2 — 8Q 2 . That is, ro cor- 
responds to the position of the unstable, circular photon orbit in the Reissner-Nordstrom 
spacetime. 
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To understand better the physics that lead to this behaviour we can use an argument 



that is due to Goebel [13]: Consider a congruence of null rays circling the black hole in the 



unstable photon orbit. To circle the black hole would require a coordinate time 

A( = 2^0 (l - ^ + V2 . (6, 

The fundamental mode frequency then follows (if the beam contains i cycles) from 

io « — . (7) 
At v ; 

In a similar way, one can infer the damping rate of the quasinormal mode from the decay 

rate of the congruence if the null orbit is slightly perturbed ]n| . From this information it 

is easy to convince oneself that the oscillation frequency of the modes should increase as Q 

increases. The results for the first few modes in Fig. [I] agree nicely with this description. 

A remarkable fact, that is notable in Fig. [TJ, is that the dipole frequencies for ^ approach 

the quadrupole frequencies for as the black hole becomes extreme. This effect was first 



noted by Onozawa et al. [[U| (for a study of the spin-3/2 case, see |]H|). They also showed 
that this surprising phenomenon arises because the corresponding two effective potentials 
are related. That is, the effective potential for i and is related to that for I + 1 and ^>2- 
Specifically, \l/i(£) corresponds to 1 — 1)- As is easily verified this is true also for the 
polar potential, but in that case the relation does not take the simple form of Eq. (C3) in 

0- 

C. Highly damped modes 



While the behaviour for the first few modes have been studied in detail before []K|-[l2 



there have been no studies of the highly damped modes for Reissner- Nordstrom black holes. 
Thus, the results for n = 5—8 in Figs. Q and [] are new. They are also truly remarkable. While 
the behaviour is rather simple for the first modes (n = — 3), the rapidly damped modes 
show many strange features. Our calculations have discovered a small zoo of seemingly 
different species. Let us discuss them separately: 



The first surprising feature occurs for n = 4 (see Figure |l|) at a charge Q ~ 0.938M: 
The ^2 mode goes through a tiny loop that closes at Q ~ 0.962M. For the following 
mode (n = 5) this loop has grown. Several similar loops are obvious in the results for 
n = 6 and ^i, and when one continues up the spectrum one finds that these loops are 
a common feature of many modes. One interesting aspect of this result is that black 
holes with different charge may share a specific mode-frequency. 

For a few modes the mode-frequencies approach the negative Im uj axis as the charge 
varies. Then our numerical methods |]il~l , |T2f become unreliable. This is the reason 
for the gaps in the data for n = 6 for and n = 7 for both ^ and \&2- If these 
modes cross the axis or not is impossible to say given the available techniques. If 
they do it would be very surprising, but given the data in our figures it does not 
seem impossible. One added peculiarity is that it seems as if some modes are "multi- 
valued" . That is, for n = 6 and ^2 we find a mode both in the upper branch close to 
the imaginary axis and in the lower branch for a small range of the black-hole charge 
(0.8985 < Q/M < 0.9105). A similar behaviour can be observed for n — 7 and 
when 0.9330 < Q/M < 0.9425. At first this result seems nonsensical and obviously 
wrong, but it is confirmed by both the continued fraction method and the numerical 
integration scheme. This could indicate that it should be taken seriously. 

Our study also sheds some light on an issue that has been debated for Schwarzschild 
black holes for some time. Is there a quasinormal mode on the imaginary axis at uoM = 
—2il Some methods indicate that there should be such a mode, and that it corresponds 
to n — 8. Other methods say that no such mode exists. To resolve this issue is difficult, 
basically because none of the proposed methods is reliable close to the axis. An 
added complication is that the case ujM = —2i is a very special one. That frequency 
corresponds to a so-called algebraically special perturbation and, as was shown by 



Chandrasekhar p3[ , one can find analytic solutions to the corresponding perturbation 



equations. Although one can convince oneself that Chandrasekhar' s special solutions 
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do not satisfy the quasinormal-mode boundary conditions, the issue is not resolved. To 
prove the existence of a quasinormal mode one must study the analytic properties close 
to the mode-frequency, and this turns out to be difficult in this specific case. Anyway, 
once the black hole acquires some small charge it is clear that a mode exists (n — 8 for 
^2 in Fig- @)- This mode approaches the suggested value ujM = —2i as Q — > 0. But 
still this does not prove the existence of a mode for the Schwarzschild black hole. It is 
plausible that the mode in Fig. ^| and its symmetric counterpart in the left half of the 
complex a;-plane will coalesce at Q — 0, and perhaps the two modes that exist for the 
charged case then cancel each other in some way. It is also worth mentioning here that 
there are algebraically special solutions also for the Reissner-Nordstrom black hole. 
For \l/ lj2 these correspond to frequencies [f23[ 



<2 = -^~ m + !){£ + 2)' 



3M =f J9M 2 + AQ 2 (£ -!)(£ + 2) 



(8) 



In general we cannot see any correlation between these frequencies and the quasinormal 
modes of the charged black hole. 

• For higher overtones than those displayed in Fig. ^ the behaviour tends to be qualita- 
tively similar to the presented ones. The mode-frequency spirals into the value for the 
extreme black hole. The spirals get tighter for more rapidly damped modes. 

How are we to understand these results? Well, at present it is very hard to make much 
sense out of the data presented in Fig. ^|. The obvious answer is that the change in the mode- 
frequencies occurs because the effective potential changes as the charge increases. For the 
slowly damped modes it is straightforward to show that the oscillation frequency depends 
on the height of the peak of the potential, while the damping rate is related to the second 
derivative of V at the peak [K| . But it is much harder to draw similar conclusions for the 
rapidly damped modes. Present methods can reliably calculate the mode-frequencies (and 
also account for the excitation of the modes in a dynamical process P0|-|2"2|]), but we have 
no clear understanding of the relation between the highly damped modes and the details of 
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the effective potential. On the other hand, we have seen that the mode-frequencies change 
dramatically as the charge increases from Q « 0.9M. If one studies the corresponding change 
in the effective potential, cf. Fig. |3], one can infer that as one approaches the extreme black 
hole i) the potential does not change much for r 3> r , ii) the peak of the potential stays 
almost constant, and iii) the main change in the potential is the fall-off rate as r» — > — oo. 
This suggests the possibility that the high overtones of the black hole are in some way related 
to the effective potential close to the horizon. This issue will have to be studied further in 
the future. 



III. CONCLUDING REMARKS 

In this paper we have presented the results of a detailed study of the quasinormal modes 
of a Reissner-Nordstrom black hole. Our work complements previous studies in several ways. 
For the first time we have considered i) the rapidly damped modes, and ii) how the mode- 
frequencies approach the anticipated value for the extremely charged black hole. We have 
given special attention to the mode-behaviour near the extreme limit. The two independent 
algorithms (from |lT| and [12[]) that we used for the calculations provide consistent results 



even for a nearly extreme black hole. These values are also in nice agreement with the results 
obtained for the extreme black hole |13 . That is, the quasinormal modes converge towards 



the values for the extreme 

We find that while the slowly damped modes can be understood (for example, by a WKB 
argument), the behaviour of the rapidly damped modes as the charge increases is harder 
to explain. In general, the modes tend to spiral into the extreme value. Moreover, there 
are cases where the quasinormal modes seem to be double-valued (in cases when the mode- 
frequency approaches the imaginary cj-axis). At one level, these are just peculiar results 
without much relevance, but they lead to several questions. For example; Will the modes 
of a Kerr black hole show similar features? This does not seem unlikely, because the Kerr 
problem is in many ways similar to the Reissner-Nordstrom one. Also in the Kerr case will 
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there be two horizons that merge as the black hole becomes extreme. Although the slowest 



damped modes of the Kerr black hole have been calculated ||15|| , the rapidly damped ones 
have not yet been considered. Since rotating black holes should be astrophysically relevant 
a study of the rapidly damped modes for Kerr could yield physically interesting results. 

The present study also leads to questions regarding the quasinormal modes themselves: 
Which features of the effective potential govern the behaviour of the rapidly damped modes? 
The present evidence (cf. Fig. seems to indicate that the fall-off rate towards the event 
horizon plays an important role. But more work is required to answer the question properly. 
At present we have to conclude that we still do not have an intuitive understanding of the 
origin of the rapidly damped black-hole modes. 

A final question regards the relevance of the algebraically special modes. In general, 
our results do not indicate any relation between the algebraically special mode (that has 
a purely imaginary frequency) and the quasinormal modes for a charged black hole. But 
it is well-known that the algebraically special mode coalesces with the ninth gravitational 
quasinormal mode in the Schwarzschild limit. Is there a profound reason for this or is it 
just a coincidence? To answer this question is difficult because all present methods for 
quasinormal modes break down in the region close to the imaginary frequency axis. On the 
other hand, it is known that the algebraically special modes of the Kerr black hole move 
away from the imaginary axis. Hence, it seems likely that a study of the Kerr problem may 
prove illuminating also for this issue. 

Several questions concerning quasinormal modes thus remain to be resolved. The present 
work has contributed a few new pieces to the puzzle, but many other pieces need to be 
added before the picture becomes completely clear. Still, the present results should help us 
get a better understanding not only of the Reissner-Nordstrom black hole, but also of the 
quasinormal- mode phenomenon in general. 
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FIGURES 




FIG. 1. The behaviour of the first six quasinormal mode frequencies (we show Im ujM as a 
function of Re tuM) for a Reissner-Nordstrom black hole as the charge is increased. The solid curves 
correspond to I = 2 and Vl/2 (which reduces to pure gravitational waves in the Schwarzschild limit) 
and the dashed curves to i = 1 and \&i (that limits to pure electromagnetic waves). To show that 
the mode-frequencies change dramatically for Q > 0.9M we have indicated the charge of the black 
hole by diamonds [at increments in Q of 0.1M]. The frequencies generally move anti-clockwise in 
the figures as the charge is increased. 
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FIG. 2. The behaviour of the rapidly damped modes (we show Im loM as a function of 
Re ujM) of a Reissner-Nordstrom black hole as the charge is increased. The right frames (labelled 
'g') correspond to I = 2 and ^2 (which reduces to pure gravitational waves in the Schwarzschild 
limit) and the left frames (labelled 'e') are for i = 1 and ^1 (that limits to pure electromagnetic 
waves). To show that the mode-frequencies change dramatically for Q > 0.9M we have indicated 
the charge of the black hole by diamonds [at increments in Q of 0.1M]. The frequencies generally 
move anti-clockwise in the figures as the charge is increased. 
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tortoise coordinate 

FIG. 3. The effective potential for £ = 2 and ^2 (that corresponds to gravitational pertur- 
bations in the Schwarzschild limit) is shown for Q/M = 0.9 — 1.0. The potential falls off slower 
towards r* = —00 as Q increases. The behaviour is similar for other values of £ and also for the 
potential that governs 
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